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ABSTRACT 



Collisions of mm-size dust aggregates play a crucial role in the early phases of planet formation. It is for example currently unclear 
whether there is a bouncing barrier where millimeter aggregates no longer grow by sticking. We developed a laboratory setup that 
allowed us to observe collisions of dust aggregates levitating at mbar pressures and elevated temperatures of 800 K. We report on 
collisions between basalt dust aggregates of from 0.3 to 5 mm in size at velocities between 0.1 and 15 cm/s. Individual grains are 
smaller than 25 jim in size. We find that for all impact energies in the studied range sticking occurs at a probability of 32.1 ± 2.5 % 
on average. In general, the sticking probability decreases with increasing impact parameter. The sticking probability increases with 
energy density (impact energy per contact area). We also observe collisions of aggregates that were formed by a previous sticking of 
two larger aggregates. Partners of these aggregates can be detached by a second collision with a probability of on average 19.8±4.0 %. 
The measured accretion efficiencies are remarkably high compared to other experimental results. We attribute this to the relatively 
large dust grains used in our experiments, which make aggregates more susceptible to restructuring and energy dissipation. Collisional 
hardening by compaction might not occur as the aggregates are already very compact with only 54 % ± 1 % porosity. The disassembly 
of previously grown aggregates in collisions might stall further aggregate growth. However, owing to the levitation technique and the 
limited data statistics, no conclusive statement about this aspect can yet be given. We find that the detachment efficiency decreases with 
increasing velocities and accretion dominates in the higher velocity range. For high accretion efficiencies, our experiments suggest 
that continued growth in the mm-range with larger constituent grains would be a viable way to produce larger aggregates, which might 
in turn form the seeds to proceed to growing planetesimals. 

Key words, planets and satellites: formation - protoplanetary disks 



1. Introduction 

Terrestrial planet formation is assumed to begin with the ag- 
gregation of micrometer-sized dust particles in protoplane- 
tary disks. This initia l growth is governed by cohesive forces 
dBlum & Wurml |2008l) . The relative velocities between parti- 
cles are determined by the different types of gas-grain coupling 
dWeidenschilling|ll977l) . At low collision velocities on the order 
of a few mm/s, fractal du st aggregates form dBlum et al.lll996t 
iDominik & Tielensll 19971) . This first phase can be characterized 
as a hit-and-stick regime, a s all collisions lead to rigid stick- 
ing at the first contact point dWurm & Blum|[l998t iBertini et al.l 
120091) . Eventually, energies become high enoug h for restructur- 
ing an d porous but non-fractal aggregates form dBlum & Wurml 
120001). Collision velocities then increase with aggregate size 
( Weidenschilling & Cuzzil["l993l) . At a certain stage, mm-sized 
compact dust agglomerate s of porosities on th e order of 60 - 
70 % might have formed dWeidling et al.ll20 09). It is currently 
unclear how planet formation proceeds from here. Depending 
on the disk model, the relative velocities rapidly increase from 
10~ 3 m/s or 10~ 2 m/s to several tens of m/s at dm to m size 
dDeschll2007l IWeidenschiiling & Cuzzilll993l) . Numerical mod- 
els cover the high speed parts and, depending on the assumed 
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parameters such as porosity, elastic propert ies or collision ve- 
lociti e s, further growth is or is not possible ( Geretshauser et all 
1201 U ISchaferetal.ll2007HWadaetaT]l2009l) . However, it an- 
pears non-trivial to reach the high speed collisions that partly 
permit growth. Recent experiments have shown that the colli- 
sion results for m m-sized dust aggregate s are more comp l ex and 
growth can occur dWeidling et al.ll201 ll) . IWeidling etall (1201 ll) 
report a sticking probability of 5.6 % at collision velocities of be- 
tween 3 cm/s and 10 cm/s that however decreases with velocity. 
Experimental studies in general have so far foun d that bounc- 
ing is a typical outcome after the sticking phase dGtittler et al.l 
l2010h . With the assumption that colliding particles bounce off 
each other, further growth would be pr evented. This assumption 
was introduced by IZsom et al.l d2010l) as the bouncing barrier. 
The underlying physics is that these collisions are more or less 
elastic and not enough energy can be dissipated by the compact 
aggregates to allow sticking. This is also the reason why growth 
is again possible at higher collision velocities as fragmentation 
leads to enough dis sipation to allow mass transfer and net growth 
dWurm et al.ll2005l) . 

One promising way to overcome the problems of collisional 
growth is the trapping of solid particles by turbulence. If the 
local particle density in the protoplanetary disk is increased 
sufficiently, gravitati onal instabilities lead to the rapid forma- 
tion of planetesimals dJohansen et a l. 2007). Particles might also 
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be co ncentrated by a streaming instability (lYoudin & Johansenl 
120071) . Current models describing this growth mechanism are 
based on the assumption that most of the solid material in proto- 
planetary disks is trapped in macroscopic (decimeter size) dust 
agglomerates. This size range therefore has to be reached by 
aggregation. If bouncing prevents growth, then these instability 
scenarios are also unable to solve the problem of planet forma- 
tion. 

In experimental studies, it has been shown that the collision 
characteristics change sig nificantly, when bodies of different 
size collide with each other. Wurm et "aT] (l200l , lTeiser & Wurml 
d2009bl) . and iTeiser et all d2011bl) demonstrated that the frag- 
mentation of the smaller body can lead to a partial mass trans- 
fer from the smaller to the larger body if the impact velocity is 
high enough for projectile fragmentation. Once the aggregates 
have reached a size of a few centimeters, the reaccr etion of dust 
by gas drag can efficiently enhance the growth rate dWurm et alJ 
I200U ITeiser & Wurml l2009at ITeiser et al.|[201 lal) . Mass trans- 
fer by (projectile) fragmentation as well as growth by reaccre- 
tion require a certain amount of aggregates to be already in the 
centimeter range. To overcome the bouncing barrier, it remains 
uncertain how at least a few particles might cross the bouncing 
b arrie r to form cm-siz e aggre g ate s . 

In IWeidling et all J201ll) . the monomer size was between 
0.5 /mi and 10 /mi (wit h 80 % between 1 /mi and 5 /mi). The 
experimental studies of iBlum et al.l d2006l) showed that the me- 
chanical properties of dust aggregates are determined by the size 
distribution of the monomers. Depending on the monomer size, 
the restructuring of dust aggregates during collisions will take 
place at varying collision velocities. Experimental evidence that 
the size distribution of the monomers influe nces the outcome 
of agg regate collisions was also presented by lLangkowski et alJ 
(2008). Although the focus of this study was on collisions be- 
tween highly porous aggregates of different sizes (mm vs. cm), 
it was clearly shown that thresholds for sticking and/or fragmen- 
tation depend on the size distribution. 

Here, we present an experimental study of the collisions be- 
tween dust aggregates with sizes between 0.3 and 5 mm at col- 
lision velocities of between 0.1 cm/s and 15 cm /s. This param- 
eter ra nge is comparable to the experiments of We idling et alJ 
( 2011b. equaling that of the bouncing barrier proposed by 



IZsom et ail ( 20101). In com p arison to previous a g gregation stud- 
ies dBlum & Wurml 120081: iGiittler et alJ l2010t Iweidling et aT] 
1201 %, the monomer size is significantly larger. Within this study, 
we use basalt with a broad size distribution of between 0. 1 /mi 
and 25 /im. As the basalt consists mostly of silicates, we re- 
gard this material as a suitable analogue for protoplanetary dust. 
In Fig. Q] we present a scanning electron microscope (SEM) 
image of the dust. The aforementioned size range for basalt 
contains the lar gest of the grains regularly found in meteorites 
dBrearlevlll9"99l) . Some authors argue that these are of nebular 
origin, t hough others also n ote that these large grains form on as- 
teroids dScott & Krolll2005]). Particles found on asteroid Itokawa 
are within the size range dNakamura et alJl201 1[) and larg e par- 
ticles were also found on comet Wild 2 dRietmeijerl l2008). With 
the exception of processed particles such as chondrules, tens 
of micrometer are at the larger end of the particle size range. 
However, this extreme might hold a key to providing the seeds 
for the growth of larger aggregates. 

2. Experimental setup 

The experimental setup is shown in Fig. [2] Within a vacuum 
chamber, a heater (up to 800 K), and an illumination source 




25 um 



Fig. 1. Scanning electron microscope image of the basalt sample 
used. The mass is dominated by monomers of size > 20//m. 



copper ring with 
LEDs (water-cooled) 




heater with 
ceramic plate 
and quartz lens 



Fig. 2. Experimental setup. A heater is placed within a vacuum 
chamber. For contrast enhancement a ceramic plate and for con- 
finement a concave quartz lens are placed on the heater. Dust 
aggregates levitate above the lens and are observed by a camera. 



(LED ring) are installed. A camera observes the collisions from 
above. To enhance the contrast, a white ceramic plate is placed 
on the heater. To confine the colliding dust aggregates to the field 
of view, a slightly concave quartz glass lens is placed on top of 
the ceramic plate. 

Basalt aggregates of millimeter size, composed of particles 
< 25 /mi, are placed loosely onto the glass. In general, fewer 
than 15 aggregates are placed simultaneously on the lens. The 
vacuum chamber is then evacuated. Experiments were carried 
out at ~ 20 mbar and ~ 800 K. In these conditions, the dust ag- 
gregates levitate. This m echanism was discovered only recently 
dKelling & Wurml 12009). While gravity is still present, there is 
no interaction (friction, sticking) with any surface and the inter- 
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X. thermal creep p 




Fig. 3. Principle of the Knudsen levitation. The pores of the ag- 
gregate act as a collection of micro-channels. Thermal creep is 
induced by the temperature difference over the aggregate. The 
pressure increases below the aggregate, and the aggregate is then 
levitated. The levitation height is limited by the gas flow to the 
sides, which r educes the overpressur e once the aggregate has 
been levitated dKelling & Wurmll2009h . 




Fig. 4. Example of levitating dust aggregates (19 mbar and 798 
K heater temperature). 



action between different aggregates can be observed. Levitation 
is induced by a pressure increase due to thermal creep through 
the pores of the aggregates. Fig.[3]illustrates the principle of the 
Knudsen levitation mechanism and Fig. [4] shows a side view of 
lev itating aggregate s . 

I Knudsen! d 19091) found that if two gas reservoirs at different 
temperatures are connected by a capillary tube whose diameter 
d t is small compared to the mean free path A of the gas molecules 
(d t <s A), the pressure increases at the warmer side. In equilib- 
rium, the pressure ratio is given by pilpx = VTT/Ti", where p2 
and Ti are the pressure and temperature in the warmer cham- 
ber and p i and T\ are t he pre ssure and temperature in the colder 
chamber. iMuntz et al.l d2002l) showed that at Kn = A ■ d^ 1 =* 1 
the pressure difference Ap = \p\ - pz\ is 



Ap 



Qt AT 

Qp T avg 



(1) 



where p avg and r avg = {T\ + 7\)/2 are the average pressure and 
temperature, Qj/Qp is the ratio of the coefficients of the thermal 
creep to back flow of the gas and Ar = \T\ - T 2 \ is the tempera- 
ture difference. 

The lateral extensions of the dust samples are much smaller 
than the glass cavity. Hence, the bottom temperature of the dust 
aggregates is determined by the heater temperature 7Y At the 
top, the aggregates cool by thermal radiation. The loss of heat is 
counterbalanced by thermal conduction through the aggregates. 
In equilibrium, the temperature T\ at the top is given by 



0-77 



7Wj 



(2) 



aggregate's vertical extension. The dust aggregates have pores, 
hence the aggregates act as a collection of micro-channels. At 
the pressures used in the experiments, the mean free path of gas 
molecules is comparable to the mean pore size of the aggregates. 

The micro-channels formed by the pores of an aggregate, 
have different temperatures at their inlet (top of the aggregate) 
and outlet (bottom of the aggregate). According to Eq.(Q]), ther- 
mal creep leads to a pressure increase Ap below the aggregate. 
If the force F^ — Ap A acting at the bottom side A of the aggre- 
gate compensates gravity Fq, the aggregate is lifted. With rising 
height, gas is released to the open sides below the aggregate. 
The levitation height is limited to some tens to a hundred of mi- 
crometers as then no additional overpressure can be established 
because the gas below the aggregate can escape to the sides. As 
the aggregates used in this experiment have typical vertical ex- 
tensions » 100/um and the levitation heights are more or less 
equal, the setup can be treated as two-dimensional. The strength 
of the levitation can be estimated as follows. We consider an 
aggregate with an area of A = 4 mm 2 , a vertical extension of 
/ = 1 mm, and a density of 1332 kg m~ 3 , which results in an ag- 
gregate mass of 5.33 mg. The pressure increas e below the aggre- 
gate c an be calculated in an analogous way to lKelling & WurrrJ 
(120091) . Using spherical, 20 pm sized particles arranged in a sim- 
ple cubic lattice, the resulting micro-channels have a minimum 
diameter of d t = 8.3 pm. With the mean free path of the gas 
molecules of A = 932 pm at a 780 K heater temperature and 
19 mbar pressure, the ratio of the coefficie nts of the therma l 
creep to back flow of the gas is Qj/Qp = 0.22 dSone et al.ll990l) . 
The top temperature of the aggregate is given by Eq.(f2]i, result- 
ing in Ti = 667 K. With 77 avg = (J x + T 2 )/2 = 724 K and 
/?avg ~ pi = 19 mbar, the induced overpressure below the aggre- 
gate (Eq. ([]])) is Ap = 0.65 mbar. The force acting on the aggre- 
gate from below is F^ - Ap A - 2.6 x 10~ 4 N. With F c = mg 
and g = 9.81 m/s 2 , the ratio of the lifting force to gravity is 



Fkn _ , 

F G 



3. Experiments 



(3) 



where cr = 5.67 x 10 8 W m 2 K 4 is the Stefan-Boltzmann 
constant, k v is the aggregate's thermal conductivity, and / is the 



To characterize the single aggregate collisions, a set of parame- 
ters of every collision is determined: the masses of the colliding 
aggregates m\ and n%2, the aggregate shapes (described by the 
circularity c), the collision velocity v, the maximum contact area 
A c between the colliding aggregates, the impact parameter I, and 
resulting parameters such as either the impact energy E or the 
impact energy density £d (defined below). 

Aggregate masses - The projected aggregate areas were mea- 
sured directly from the two-dimensional (2D) video recordings. 
Based on these measurements we model the aggregate volume 
by assuming a spherical hemisphere (steep rise to the edge in 
Fig.|5]l. We compared the resulting vertical structure of the par- 
ticles with side view images from the aggregates. The model 
reproduces the shape and hence the mass to within an accuracy 
of about 50 %. 

Circularity - The aggregates used in the experiments are far 
from being spherical. The deviation of the aggregate shape from 
the spherical shape is expressed by the circularity c which, is 
defined as 

(4) 

where A agg is the projected aggregate area, A max = n r 2 ^ is the 
area of a circle and, r max is the maximum distance between the 
center of mass and the edges of the aggregate (Fig. [6j». 
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Fig. 5. Example of an aggregate volume approximated by a 
hemisphere (arbitrary units). 
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Fig. 7. Example approach of two colliding aggregates with a col- 
lision velocity at t — Oms of v ~ 1.8 cms -1 . The movement is 
approximated by a second order polynomial (solid line). In com- 
parison, a linear fit to the first four data points is shown to visu- 
alize the decelerating component of the sideflow of gas, which 
is induced by the Knudsen effect. 



t t 



Fig. 8. Example of a collision between two aggregates with a 
time step of ~ 8 ms from the left to the right. The red (gray) line 
corresponds to the maximum contact line l c . 



'max 



Fig. 6. The circularity c of an aggregate is defined by the relation 
of the aggregates area A agg (dark-gray) and the maximal area 
A max (light-gray) of a circle defined through r max , which is the 
maximum distance between an aggregate pixel and the center of 
mass (black cross). 



Impact velocity - The impact velocity v is calculated from 
the distances between the aggregate centers of mass. The centers 
of mass of the aggregates are calculated directly from the dis- 
crete pixels of the projected aggregate areas via r com = r,-« -1 , 
while r, are the respective coordinates of the aggregate pixels 
and n is the total number of aggregate pixels. Owing to the gas 
outflow on the sides below the aggregates, the approach of two 
aggregates is slightly decelerated. We used a second order poly- 
nomial to fit the approach of the aggregates. An example of this 
fit can be seen in Fig. [7] 

The collision velocity depends on - besides the accuracy of 
the fitting parameters - the exact point in time when the two ag- 
gregates touch. The collision time is defined through the contact 
of the two aggregates within the 2D image sequence with an un- 
certainty of one image or 2.5 ms. For the slowest collisions (1 
cm/s) with a low total mass, this implies a maximum uncertainty 



of 6.5 %. For fast collisions with high total masses, the decel- 
eration is less intense and the velocities are more accurate with 
uncertainties of up to 3 %. 

Maximum contact area - The contact area is important in 
particle collisions as it determines the number of sticking con- 
nections and the way in which the force is distributed among the 
grains. We consider the maximum contact area A c of the aggre- 
gates visible in the 2D images. If the maximum contact line l c 
(see Fig. [8} during a collision were equivalent to the diameter 
of a circle, the circular contact area would be nil ■ 1 /4- The ag- 
gregate's bottom side is flat because it was lifted from a plane 
surface. Hence, the assumption of a full circle as a contact area 
leads to an overestimate of the contact area by a factor of two 
giving n% ■ 1/8. As the lateral extension of the aggregates is 
usually larger than the height, we assume a more ellipsoid shape 
of the contact area, which adds a small factor of 0.8 to the max- 
imum contact area, a conservative improvement but otherwise 
arbitrary. Finally, the maximum contact area is 



10 ' 



(5) 



Impact parameter - We define an impact parameter / to de- 
scribe the difference between central and oblique collisions as 



/ = 



r\ + r 2 



(6) 



This is a dimensionless value between (central collision) and 
1 (contact only on the outer edges of both aggregates), where 
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projected trajectory of ,*** 
the smaller aggregate..** 



7H 




/ offset trajectory 
/ through the center 
of mass of the larger 
aggregate 



Fig. 9. Visualization of the components used for the impact pa- 
rameter calculation, where r\ , ro and d are determined on the last 
frame before the aggregates get in touch. The dotted line repre- 
sents the momentary projected trajectory of the smaller aggre- 
gate relative to the center of mass of the larger aggregate (white 
cross). The offset trajectory through the center of mass of the 
larger aggregate is marked with a dashed line. The distance d 
between the center of mass of the larger aggregate and the mo- 
mentary projected trajectory of the smaller aggregate is visual- 
ized with a thin line. 



d is the distance between the projected trajectory of the smaller 
aggregate and the center of mass of the larger aggregate (Fig. 
[9]). The dimensions r\ and are defined to be perpendicular to 
the impact direction one frame before the aggregates come into 
contact. 

Impact energy - The impact energy E is calculated using the 
impact velocity v and the calculated masses m\ and ma via 



1 mi mi 

E = l -^—v 2 . 

2 m\ + m.2 



(7) 



4. Results and discussion 

After the data reduction, the following parameters of a collision 
are given: the aggregate masses m\ and ni2, the circularities c\ 
and C2 of the aggregates, the impact velocity v, the impact energy 
E, the maximum contact area A c , and the impact parameter /. In 
the following, we divide the collisions into two types: 

Normal collisions - These are collisions between the aggre- 
gates that were placed initially on the heater (a total of 239). 

Pre-accreted collisions - These are collisions between aggre- 
gates that were formed before by a sticking collision (a total 
of 101). 

The latter collisions are clearly visible as two connected 
aggregates. Most basic models of dust aggregate collisions in a 
protoplanetary disk provide a collision velocity depending on 
the two particle masses (and sizes). In Fig. [10] the outcomes 
of the individual normal collisions and in Fig. [TT] those of the 
pre-accreted aggregates are plotted depending on the collision 
velocity v and the sum of the mass of the two aggregates (green: 
sticking; yellow: bouncing; red: detachment). 

The sum of the aggregate masses range between 10~' mg 
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Fig. 10. Outcome of the individual normal collisions: green in- 
dicates sticking particles and yellow those bouncing. 
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and 20 mg and their velocities between 0.3cm/s and 15cm/s. 
For all 340 collisions, sticking occurs with a probability of 
32.1 ±2.5 %. The sticking probability for normal collisions is on 

average 35.1 +3.0 % and for pre-accreted collisions 24. 8 + 5. 3 %. p(v) = a v (l - e~^ vV ) 



Fig. 11. Outcome of the individual pre-accreted collisions: green 
sticking particles, yellow those bouncing, and red detachment. 



The probability that a collision between two pre-accreted aggre- 
gates leads to detachment is 19.8 + 4.0 %. In the following, we 
evaluate how the sticking probability depends on the different 
parameters mentioned above. In general, we give three proba- 
bilities, one probability for sticking p s , one for bouncing p\> and 
one for detachment p&. The number of data points in a bin is 
typically 25, depending slightly on the chosen data set. The er- 
ror bars in the probabilities mark the standard deviation in the 
mean. Error bars of certain quantities also include measurement 
errors. At the original spatial resolution of ~ 50 /^m/pixel, we 
do not observe splits of dust aggregates in collisions but at high 
spatial resolution (~ 7.5 /^m/pixel) we see small fragments or 
small-scale mass transfer from one aggregate to another (Fig. 
[T2l . Partial compaction of aggregates was also observed in the 
high-resolution configuration. 

The velocity dependence of the sticking probability of the 
239 normal and the 101 pre-accreted aggregate collisions are de- 
picted in Figs[l3]and[T4] There is a clear tendency for the stick- 
ing probability to increase with collision velocity. We used the 
following analytic expression to fit the data for normal collisions 



(8) 
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1 mm 




Fig. 12. High spatial-resolution record of a collision with a mass 
transfer at as 5.2 cm/s. A small fragment from the right aggregate 
is attached to the left aggregate (circle) after the collision. 




4 5 6 

velocity [cm/s] 

Fig. 13. Velocity dependence of the sticking probability of nor- 
mal collisions. The dotted line corresponds to p s (v) = 0.47 ■ (1 - 
""")■ 



-0.32 i 



with a v = 0.47 ± 0.07 and p v = (0.32 ± 0.10) s ctrT 1 . The reason 
behind this choice is as follows: for small velocities, we expect 
bouncing as no restructuring occurs and the collisions are elas- 
tic. For large velocities, the sticking probability has to level off 
at a value smaller than one by definition. Equation © is a simple 
function that fulfills these requirements. For pre-accreted aggre- 
gates, detachment also occurs. At higher velocity, the detach- 
ment probability is significant lower however than the sticking 
probability. 

The overall (normal and pre-accreted) sticking probability 
dependence on the collision energy is shown in Fig. Q3] There 
is only a slight trend that higher energy leads to an increase in 
the sticking probability, but for irregular particles the contact ar- 
eas vary and play an important role in the energy distribution and 
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sticking 
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velocity [cm/s] 



10 



12 



Fig. 14. Velocity dependence of the sticking probability of pre- 
accreted collisions (dark-gray: accretion; gray: bouncing; light- 
gray: detachment). The colored bars represent the probabilities 
and were placed at the average velocity of the respective bins, 
which are indicated by the vertical lines. 




o.i 

energy [nJ] 

Fig. 15. Impact energy dependence of the sticking probability for 
normal and pre-accreted collisions. 

dissipation. We therefore define an energy density as - E/A c , 
where E is the impact energy and A c is the maximum contact 
area. The larger his value, the more the restructuring that can oc- 
cur to dissipate the energy that might be beneficial for sticking. 
Figure Q~6] indeed shows a clear increase in the sticking proba- 
bility (normal and pre-accreted) with energy density in the in- 
vestigated regime. We used the same function (Eq. [9]l to fit the 
data 



p(£d) = OE d (l-e^ £j ) 



(9) 



witho-E, =0.40 + 0.02 <md/3 Ed = (13.4 + 2.6) mm 2 nJ" 1 . 

Collisions can also be characterized by the mass ratio of the 
impacting particles. The sticking probability for normal colli- 
sions is depicted in Fig. [17] Within the studied interval from 
equal masses to a ratio of 100, a significant dependence is not 
visible. 

To see a potential shape effect, we show the sticking effi- 
ciency on the average circularity c avg = (c\ + c£)l2 for normal 
collisions in Fig. [18] There is no trend in the data but instead a 
rather large variation. 
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energy density [nJ/mm 



Fig. 16. Impact energy density dependence of the sticking prob- 
ability Ed of all collisions, where the appropriate function is 
/? s OE d ) = 0.40-(l- e - 13 - 4 H 




Fig. 17. Sticking probability versus mass ratio for normal colli- 
sions. 

The last dependence considered is on the impact parame- 
ter I, which in general is known to be important in collisions. 
The probabilities are plotted versus the impact parameter for 
both normal collisions in Fig. [19] and pre-aggregated collisions 
in Fig. [20j where a clear dependence is visible. For normal colli- 
sions (Fig. |T9T>, the sticking probability first decreases to a min- 
imum at I as 0.2. This could be because head-on collisions can 
easily lead to deformation but with increasing impact parame- 
ter these deformations were diminished. The sticking probability 
then increases to a maximum at / « 0.3. We attribute this to the 
possibility of dissipating energy by inducing particle rotation. 
The sticking probability then decreases toward larger impact pa- 
rameters because the normal velocity component gets smaller 
with less energy dissipation, while the tangential force at the 
same time is eventually sufficient to allow sliding with some 
friction but inertia to dominate. A simple analytic expression to 
quantify the decreasing part is 



Ps(I) = -J + 



1 



I + T 



(10) 



with y = 0.38 ± 0.02 and t = 0.84 ± 0.04. 

For pre-accreted collisions, the behavior is quite similar (see 
Fig. [20b . First the sticking probability increases to a maximum 




0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 
average circularity 

Fig. 18. Dependence of the sticking probability of normal colli- 
sions on the average circularity c avg . 




0.4 0.6 
impact parameter 

Fig. 19. Sticking probability versus impact parameter for normal 
collisions. The dotted line is a fit for / > 0.2 with p s (I) = -0.38 + 
(7 + 0.83)- 1 . 



at I * 0.3 and then decreases toward larger impact parameters. 
The detachment probabilities for pre-accreted collisions seem to 
be slightly larger at larger impact parameters. This would be in 
agreement with inducing rotation and tangential forces between 
sticking aggregates that would more easily break the contact. 
The effect is not strong and an increase is statistically not ob- 
served at a significant level. 

At very low velocities, repulsive gas drag might lead to a 
sticking probability that is systematically lower. In addition a 
number of detachments are likely supported by the repulsive gas 
drag as we also observed detachment without any further col- 
lision. We did not record a side view with every collision and 
aggregate and at the highest impact energies in particular aggre- 
gates might slide over each other, which would lead to sticking 
where free collisions would not necessarily result in sticking. 
Collisions for which this bias was obvious in the videos were not 
taken into account here but we cannot exclude these effects on 
a non-resolvable scale. This would lead to a systematic increase 
in the sticking efficiency. On the other hand, that the sticking ef- 
ficiency depends on the impact parameter also suggests that the 
sticking enhancements due to the two-dimensional setup do not 
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Fig. 20. Sticking probability versus impact parameter for pre- 
accreted aggregates (dark-gray: accretion; gray: bouncing; light- 
gray: detachment). 



dominate our results as other dependences should not then be so 
pronounced and agree with general expectations. 

To summarize our results, our experiments show that colli- 
sional outcomes (sticking probabilities) in the studied parameter 
range depend significantly on only three parameters: impact pa- 
rameter, collision energy density, and collision velocity. 



5. Conclusion 

There is currently a debate about whether the growth of dust 
would stal l at mm-size, whic h has been termed the bouncing 
barrier by IZsom et a D d2o!oh . This barrier would arise if par- 
ticles entered a regime of equally sized mm particles colliding 
only among themselves at velocities < 1 m/s, where only bounc- 
ing was the result of a collisi on. The idea is already present in 
early laboratory experiments dBlum & Muenchlll993l) . Rebound 
is also seen in n- particle simulatio ns depending on the coor- 
dination number (Wad a et al. I limb. Numerical modeling, in- 
cluding all different kinds of collisions then result in a bounc- 
ing barrier dZsom et al.ll2010h . However, the most pronounced 
input is that only rebound was observed in collisions of mm- 
size aggregates in the relevant domain of low c ollision velocities 
(iHeifielmann et alteOlCArWeidling et al.l20TTb . The experiments 
that are important used dust with grain sizes of about l/im. This 
is reasonable as for instance matrix ma terial in many p ristine me- 
teorites contain particles of this size dBrearlevHl999h . It is also 
reasonable as dust aggregates of this size are subject to strong 
cohesive forces that might be beneficial to growth by means of 
sticking. However, for compact mm-aggregates the strong stick- 
ing might be counterproductive for further growth. Aggregates 
are so strong in slow collisions that no energy dissipation by re- 
structuring can occur. The collisions then get mostly elastic and 
rebound is the outcome. 

That rearrangement of grains is beneficial to growth has e.g. 
been sho wn by impacts of c m-size aggreg ates at velocities o f up 
to 2 m/s jBeitz et alj|201 lb. 1.5 - 6 m/s dKofhe et al.ll2010h , or 
even 60 m/s ( Teiser & Wurmll2009bl) . These experiments show 
that as soon as the destruction of a projectile is possible, some 
mass might be transferred to a somewhat larger target, leading 
to net growth. 

Slow collision experiments with lar ger individual grains 
were carried out by iHartmannl dl978l) and lColwell eTail d2008b 



showing that larger solid projectiles might stick to the dust tar- 
get in the sub m/s range. Obviously going to a more granular 
medium with less cohesive forces might help growt h at least at 
certain stag es. In addition, the recent experiments of iBeitz et alJ 
dsubmittedl) on mm-size particles with dust rims assume that 
including larger particles helps us to dissipate energy and in- 
crease sticking velocities. In principle, varying the monomer 
sizes shifts the velocity ranges where certain collisional regimes 
such as sticking, bouncing, and fragmentation dominate. 

The grains used in our experiments are up to 25 /im in 
size. With a 54% ±1% porosity, the aggregates are already 
rather dense and unlikely to be compacted much further in low 
speed collisions. Though aggregates consisting of these parti- 
cles are fragile, they still experience sufficient sticking to form 
aggregates. What they provide by their low sticking forces is 
a mean energy dissipation in the slow collisions that leads to 
the large observed sticking probabilities. Therefore, if present, 
larger grains will aid collisional growth at low velocities, which 
might bridge growth to particles of cm-size or larger. 

Caveats - (1) If the experimental results are to be applied to 
planetesimal formation, particle size is very important. It is cur- 
rently still an unknown property of particles in protoplanetary 
disks but the basalt particles used are likely among the largest 
of possible grain sizes (excluding chondrules). Nevertheless, if 
only a few aggregates consisting of large grains formed in the 
disk just by statistical reasoning, they might have acted as seeds 
for further growth but this has to be studied in further experi- 
ments. 

(2) More fundamentally, one has to be aware that our experi- 
ments were confined to two dimensions and the aggregates were 
supported against gravity. However, the observed trends agree 
with expectations, e.g. the dependence on the 2D impact param- 
eter would also be expected in a three-dimensional case. This 
indicates that the high sticking probabilities are no artifacts of 
the experiment in general. 

We were unable to determine whether accretion predomi- 
nates over detachment based on the experiment so far and the 
bouncing barrier might still exist for most of the collisions. 
Nevertheless, sticking can be efficient and allows the formation 
of larger aggregates. By simple statistical reasoning, as not ev- 
ery aggregate becomes detached by the subsequent collisions, 
the growth of larger p articles seams feasible. The model of 
IWindmark et al.l d2012b requires some seeds but not too many, a 
situation that the collisions of aggregates with larger constituents 
might help to improve. 
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